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We calculate the photonic band gap of triply periodic bicontinuous cubic structures and of tubular 
structures constructed from the skeletal graphs of triply periodic minimal surfaces. The effect of the 
symmetry and topology of the periodic dielectric structures on the existence and the characteristics 
of the gaps is discussed. We find that the C(l2-Y**) structure with laid symmetry, a symmetry 
which is often seen in experimentally realized bicontinuous structures, has a photonic band gap with 
interesting characteristics. For a dielectric contrast of 11.9 the largest gap is approximately 20% for 
a volume fraction of the high dielectric material of 25%. The midgap frequency is a factor of 1.5 
higher than the one for the (tubular) D and G structures. 
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I. INTRODUCTION 



Photonic crystals are periodic dielectric composite structures that forbid propagation of electromagnetic waves, in 
any direction and for any polarization, for a certain frequency range, called the photonic band gap (PBG). PBGs are 
characterized by two quantities: the width of the bandgap and the midgap frequency. The characteristics of PBGs in 
periodic dielectric structures depend on the dielectric contrast between the composites, the symmetry and topology 
of the structure and the filling factor, the ratio between the volume occupied by each dielectric with respect to the 
total volume of the composite. 

Many applications of photonic crystals are in the visible (400-700 nm) or near-infrared wavelength (700-1300 run) 
range (3, El- El El ■ Finding three-dimensional (3D) photonic structures with a large absolute band gap in the 
visual or near-infrared frequency range and suitable for large-scale production is therefore highly desirable. Because 
of the sub-micrometer resolution required in the production technology, manufacturing these 3D photonic crystals is 
a great challenge. Techniques involved in the manufacturing process of the PBG materials include microfabrication, 
microlithography, self-assembly methods, weaving of dielectric fibers and the glancing angle deposition technique 0. 

In the two last decades several 3D photonic crystal architectures have been proposed and/or produced. Amongst 
them are the criss-crossing pore structures (inverse diamond lattice of overlapping air spheres in a high dielectric 
background) 0, Q , the woodpile structure 0, > inverse opal structures 0, 03 Ej| , circular spiral structures ITU , 
structures from woven dielectric fibers ^(|, s quar e spiral structures [l7L[T^ |. self-assemblies of metal nanospheres |l9l 
I20I ] and bicontinuous cubic microstructures |2ll l22j . In self-assembling systems many cubic phases with various 
symmetries and topologies have been discovered, see for example Refs. [S IS IS |2(j. However, the list of the 
experimentally identified cubic phases may be incomplete because the analysis of the scattering patterns of these 
complex structures can be very difficult [271 |28j . In self-assembling systems the gyroid structure with laid symmetry 
is observed quite often. Mesoporous solids (silicas) exhibiting these gyroid morphologies have also been obtained [2j| 
IS El- It was suggested that these structures have potential as PBG materi als |30ll . Recently, however, it has been 
shown theoretically that these double gyroid structures do not exhibit a PBG |22|. The more simple cubic structures 
such as the P (Pm3m), D (Fd3m) and G(/4i32) structures can develop PBGsjSIS]- 

The purpose of this paper is to calculate the band gaps for a large variety of triply periodic bicontinuous cubic 
structures: The P, D, G, C(P), C(D), Y, C(Y), ± Y, C( ± Y), C(Y**)=C(G), S, C(S), I-WP, F-RD, I 2 -Y**, C(I 2 -Y**), 
and L |82t |33| structures and the P, D and G tubular structures. Thereby we focus on the influence of the symmetry 
and topology of the dielectric structure on the appearance of a PBG. 

The paper is organized as follows. Section 2 describes how we generate the different periodic dielectric structures 
and how we calculate the PBGs. In sections 3 and 4 we discuss the results. Section 5 contains our conclusions. 
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II. THEORY 

In mesoscopic self-assembled bicontinuous cubic structures, interfaces separate adjacent regions of different com- 
position. These interfaces are triply periodic, i.e. they have a 3D Bravais lattice, and are often called intermaterial 
dividing surfaces (IMDS). IMDS can be approximated by constant mean curvature surfaces which can be modelled by 
level surfaces |32| . Single level surfaces, dividing space into two infinite, connected but disjunct regions, are defined 

by mm 

f(x, y,z)=J2 \F(hkl)\ cos(hX + kY + lZ- am) = t, (1) 

hkl 

where X — 2-nx/a, Y — 2ny/a, Z — 2nz/a, (x,y, z) are the positions of atoms in the crystal structure and a denotes 
the length of the crystallographic cell, the smallest cube generating space by the lattice a hki and \F(hkl)\ 

denote the phase angle and the structure factor amplitude, reflecting the symmetry of the structure, respectively. The 
parameter t determines the volume fraction of the two regions, or labyrinths. For t = 0, Eq. (1) defines the nodal 
surfaces |33|| . They are used as approximations to the triply periodic minimal surfaces, surfaces for which the mean 
curvature is zero at every point. Note that t = does not always divide space in two regions of equal volume. The 
labyrinth volume of the single level surfaces decreases monotonically with increasing absolute value of t [35|. At t\ 
and i 2 , where t\ < < £ 2 , the surfaces "pinch-off", i.e. become disconnected. They do no longer subdivide space 
into two continuous subvolumes but reduce to a lattice of closed-packed units with a given symmetry [^E ^ . The 
pinching off behavior, and hence the change in topology can be controlled by considering the skeletal graphs of triply 
periodic surfaces generated using level set methods |3jj . After the surfaces pinched-off , they disappear completely for 
some t < t\ or t > t%. 

The two regions defined by the single level surfaces are not simply connected. They interpenetrate each other 
in a complicated way. The topological complexity of the level surfaces can be characterized by the genus g [38|. 
The genus counts the number of handles, i.e. the number of holes in the closed surface. A finite surface of genus 
g is the topological equivalent of a sphere with g handles. The genus of a sphere is zero, whereas the genus of a 
torus with one hole is one. In what follows we give the genus for the crystallographic cell. The genus is computed 
using the integral-geometry morphological image analysis method |ssj . The two labyrinths may differ in shape or 
there may exist symmetry operations mapping one labyrinth onto the other. For the case of the PBG structures we 
consider the two regions to consist of material with a different dielectric constant. This makes the two labyrinths 
always distinguishable, even if they have the same shape. The symmetry properties of the level surface structures 
are therefore described by the crystallographic space group that does not include symmetry operations which would 
interchange the different dielectric regions. 

The PBG properties of dielectric structures with IMDS that can be modelled by P, D, G minimal 0, [SJ and 
level surfaces [23] have been studied previously for a dielectric contrast of 13. In this paper we report on band gap 
calculations for periodic dielectric structures defined by the single level surfaces P, D, G, C(P), C(D), Y, C(Y), ± Y, 
C( ± Y), C(Y**)=C(G), S, OS), I-WP, F-RD, I 2 -Y**, C(I 2 -Y**), and L [HIH and b y the skeletal graphs of the P, 
D and G minimal surfaces [33 • Since the I-WP, F-RD, I 2 -Y**, C(I 2 -Y**) and L surfaces for t = do not divide space 
in two regions of equal volume, we investigate both the direct and the inverse dielectric structures defined by these 
surfaces. The cubic periodic dielectric structures consist of 4 x 4 x 4 unit cells of length a = 3A, where A is the unit 
of length, and are characterized by magnetic permeability fj,(x,y,z) — 1 and dielectric constant 

six v z) - { £s ' if f( x ^ z ) <t ( 2 ) 

where £& and e s denote the dielectric constant of the background and structure, respectively. We choose the dielectric 
constant of the structure material to be that of Si (e = 11.9) jl^ and the one of the background material to be that 
of air (e = 1). The volume fraction or filling ratio of the dielectric material forming the background and structure 
are denoted by <f>b and <f> s , respectively. Varying the parameter t changes <fi s and fa and hence the composition of the 
material. The inverse structures are defined by 

•<-*•>- {Z: " %z$V*- (3) 

For the various periodic dielectric structures we obtain the band gap maps, showing the band gap width as a 
function of the volume fraction cf> s , from the computation of the density of states (DOS), or eigenvalue distribution. 
The DOS is calculated using an unconditionally stable finite difference time domain (FDTD) method to solve the 
time-dependent Maxwell equations 0, . The algorithm used is the T4S2 algorithm and the algorithm-specific 
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parameters are r = At = 0.075A/c, 6 = 0.1A, M = 8192 (see Refs. Tne width 01 the bandgap is defined 

as the width Alu of the frequency range of the bandgap divided by the midgap frequency loq and is expressed in 
percentages. Angular frequencies are given in units of 2nc/a where c is the speed of light in vacuum. 

III. LEVEL SURFACE STRUCTURES: RESULTS 

In Tables I-IV we give the equations for the direct and inverse level surface structures together with 3D renderings 
of 2x2x2 unit cubes of the structures for various values of t. Note that the equations for the direct and inverse 
structures are the same. The structures become different through Eq.(2) and Eq.(3) for the direct and inverse 
structure, repectively. The interval of t values is chosen such that the endpoints approximate the t values for pinch-off 
of the surface and such that the range of t values for which the structures show a PBG is covered. As can be seen 
from Tables I-IV the length of these intervals can be quite different for the various level surface structures. The 
corresponding volume fractions 4> s are also given. The structures are listed in the order of increasing index of their 
space group |42| (see also Table V) and not in order of appearance in the text. 

We first consider the P, D and G level surface structures. For the P surface a family of connected level surfaces 
exists for —1.45 < t < 1.45, corresponding to (f> s values between 0.14 and 0.86, where we calculated the volume 
fractions 4> s by means of the morphological image analysis method |38l ]. The D level surface family exists from a 
volume fraction <f> s of 0.16 to 0.84 corresponding to t values between and \ti\ — = 1.00. The G level surfaces 
pinch-off at \ti\ = = 1.41 which corresponds to values of <fi s of 0.02 and 0.98. 

Figd shows the band gap maps for the P, D and G level surface structures. Gaps for a broad r ang e of volume 
fractions are observed for the P, D and G structures, in agreement with the results presented in Ref. |22j. For the P 
structure a gap is found for volume fractions <j> s in the range 12-33 %, corresponding to t values between -1.50 and 
-0.60. Note that for —1.50 < t < —1.45 the structure is disconnected and consists of "spheres" on a simple cubic 
lattice. The largest gap Aui/uiq = 6.8% is found for a volume fraction of approximately 24%. The corresponding 
midgap frequency ujq — 0.442. A PBG is absent for <j) s = 50% (t = 0), in agreement with Refs. [2lll22]|. The band 
gap maps for D and G level surface structures show that compared to the P structure, the D and G structures have 
gaps for a much broader range of volume fractions. The maximum value for the D and G level structures is reached 
for volume fractions of 0.25 (t = —0.60) and 0.21 (t = —0.90), respectively. For the D (G) structure Aoj/ojq — 20.3% 
(21.7%) and loq — 0.541 (0.474). Comparison of the band gap maps in Fig.l with the ones in Ref. [23 shows that 
the PBGs in Fig.l close at slightly lower filling fractions. This may be due to the fact that we have chosen a lower 
dielectric contrast between the media. In general, lowering the dielectric contrast decreases the band gap width. In 
contrast to the numerical method used in Ref. |22j our method puts accurate bounds on the gap, even at low volume 
fractions. Hence the band gap maps in Fig|T]for the P and D level surface structures close, in contrast to the ones 
shown in Ref. |22| . The maximum value for the bandgap widths, the midgap frequency and the corresponding filling 
fractions are summarized in Table V. 

As mentioned above the characteristics of a band gap not only depend on the dielectric contrast and the volume 
fractions of the two media but also on the symmetry and the topology of the structure. The P, D and G level surfaces 
have Pm3m, Fd3m and J4i3 2 sy mmetry, respectively [UES- The P surface has genus 3, the D surface has g = 9 
and the G surface has g = 5 [38|. Thus, it is of interest to repeat the band gap analysis for level surface structures 
having the same symmetry but with different topologies. 

We first consider the space group Pm3m jSJ of the P surface to which also the C(P) and TWP surfaces belong [3^, 
l43l ]. The C(P) level surfaces with g = 9 pinch-off at \ti\ = |ta| = 0.67 which corresponds to values of (j) s of 0.28 
and 0.72. These volume fractions are relatively high compared to the volume fraction at which the largest PBG 
for the P structure occurs. For the I-WP surface, which has genus 7, a family of connected level surfaces exists for 
—3.00 < t < 2.98, corresponding to (j) s values between 0.10 and 0.99. Band gap calculations for both structures show 
no gap at any volume fraction for t\ < t < ti- Since for the I-WP surface, t = does not divide space in regions of 
equal volume, we also study the inverse TWP structure. Also for this structure we find no PBG. Hence, for the level 
surface structures with Pm3m symmetry only the structure with the lowest genus, and hence the simplest topology, 
has a PBG. 

A level surface which has the same symmetry as the D surface is the C(D) surface [33Ll43| . The C(D) level surface 
family only exists for relatively high volume fractions: 0.35 < (j> s < 0.58 corresponding to —0.28 < t < 0.26. The 
C(D) surface has genus 121. For the C(D) surface, two different values for g for the crystallographic cell can be found 
in the literature: g = 145 as found by Fisher et al. ^3 f° r example, and g = 121 as found by Garstecki et al. [44| . 
This difference might be caused by the fact that Fisher et al. have worked with minimal surfaces, while Garstecki et 
al. have used the nodal approximations, as we do in this work. The same reasoning can be applied to explain the 
fact that the C(D) level surface does not divide space in regions of equal volume for t — 0, while the minimal C(D) 
surface does [43j- Our calculations indicate that the C(D) structure shows no PBG. 
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The C(Y**)=C(G) surface, also called the complementary gyroid [3j|, has the same symmetry as the G surface 
and has g — 29. The C(Y**)=C(G) level surface family with g — 29 exists from a volume fraction <fi s of 0.29 to 0.71 
corresponding to t values between 0.00 and \t\ = 2.14. For larger \t\ values topologically new intermediate surfaces are 
formed before the pinching-off of the surface. For \t\ values between 2.15 and 3.32 the surface changes into a surface 
with genus 5, the same value as for the G surface. For \t\ values between 3.33 and 3.40 the surface changes again into 
a surface with genus 29. The surface finally pinches off at |ii| = \t2\ = 3.41, corresponding to volume fractions of 0.16 
and 0.84. As can be seen from FigEl the C(Y**) structure exhibits a gap for volume fractions 4> s in the range 11-38 
%, corresponding to t values between -4.00 and -1.25. However, the largest part of the band gap map, namely for 
0.17 < <j>s < 0.29, shows gaps for surfaces which have g = 5 and not g = 29. The maximum bandgap width of 11.9% 
is reached for cj) s — 0.18 (t = —3.25). The corresponding midgap frequency loq = 0.514. For this volume fraction the 
structure is broken apart and consists of 20 separate units per unit cell. 

For the space groups Pm3m and Fdim only the level surface structures with the lowest genus, namely the P and D 
structures, respectively, have a PBG. For the I4i32 space group, the G structure with g = 5 and the C(Y**) structure 
with g = 29 have a PBG. The gap of the G structure, however, is much larger than the one of the C(Y**) structure 
and is thus far more interesting. Moreover the largest part of the band gap map of the C(Y**) structure is not formed 
by surfaces which have g = 29 but which have g = 5, the same genus as the gyroid. Hence, for the space groups 
Pm3m, Fd3m and I4i32, the cubic structures with the simplest topology (lowest genus) have the largest PBG. In 
what follows we investigate whether this is a generic property or not. We study the band gap maps of the bicontinuous 
cubic structures in the order of increasing index of their cubic space groups [43] (see Table V). 

Two surfaces having the symmetry Pa3 but having different genus are the ± Y and the C( ± Y) surfaces. For the 
± Y structures with g — 21, which exist for 0.35 < <f) s < 0.65 (|ii| = fa] = 0.37), no gap is found for the range of 
parameters analyzed. However for the C( ± Y) structures we find a band gap. The variation of the band gap with 
volume fraction for the C( ± Y) structures with g — 13 is shown in Fig|21 The existence region of the C( ± Y) level 
surface structures is 0.07 < 4> s < 0.93, corresponding to —1.65 < t < 1.65. A PBG is found for volume fractions 
between 0.08 and 0.37 (—1.60 < t < —0.50) and is largest for a volume fraction of 17% (t = —1.30). The largest band 
gap (8.1%) is not so large as the one for the D and G surfaces but the midgap frequency uio — 0.786 is noticeable 
larger. This can be advantageous for the production of PBG crystals in the shorter wave length regimes. Also for 
this space group only the dielectric structure with the lowest genus shows a gap. 

The next cubic space group we consider is the 7a3 group. The C(S) surface belongs to this group. For the nodal 
approximation of the C(S) surface the morphological image analysis method [38| yields g = 65. This value for the 
genus is much higher than g = 17 reported by Fisher et al. for the minimal C(S) surface For the C(S) level 

surface family connected level surfaces exist for —1.41 < t < 1.41, corresponding to <fi s values between 0.23 and 0.77. 
Our calculations show that the C(S) structure has no PBG, which may be not so surprising due to the high genus 
and hence complex topology of the structure. 

Two other level surfaces belonging to a cubic symmetry group based on a simple cubic Bravais lattice are the Y 
and C(Y) surfaces with P4a32 symmetry. Both level surfaces have g = 13. Fisher et al. found g = 17 for the Y 
surface [43] . This difference might again be caused by the fact that Fisher et al. have worked with minimal surfaces, 
while we use the nodal approximations. The Y level surface structure, which exists for volume fractions between 
0.27 and 0.73 (—0.44 < t < 0.44), has no PBG. For the C(Y) surface, a family of connected level surfaces exists for 
— 1.84 < t < 1.84, corresponding to volume fractions in the range 8-92%. The band gap map for the C(Y) structure is 
shown in FigEl The C(Y) structure exhibits a gap for volume fractions ranging from 0.16 to 0.35 (—1.50 < t < —0.70). 
The maximum value Aw/wo = 6.2% (uio = 0.630) is reached for a volume fraction of approximately 22% (t = —1.25), 
as can be seen from Table V. The Y and C(Y) level surface structures have the same symmetry, the same genus and 
the same dielectric contrast but only the C(Y) structures have a PBG. The Y structures, however, only exist for a 
relatively high fraction of the volume occupied by the high dielectric material, which in general is a disadvantage for 
the appearance of PBGs. 

A single level surface belonging to a cubic space group based on a body-centered cubic Bravais lattice is the S 
surface. The S surface has /43<i symmetry and has genus 21. The family of connected S level surfaces exists for 
—0.74 < t < 0.74, corresponding to <fi s values between 0.12 and 0.88. Our calculations show that the S structure has 
no PBG for the parameters studied. 

The band gap maps of the F-RD and the inverse F-RD structures are depicted in Fig|3| The F-RD surface has 
FmZm symmetry and genus 21. The cubic space group Fmim is a supergroup of the group Pm3m, the symmetry 
group of the P surface [33 . |43| . The groups differ only by a face-centering translation. The F-RD surfaces pinch off at 
ti = —1.12 and t^ = 0.99 which corresponds to values of <f> s of 0.21 and 0.67. The results show that the F-RD structure 
exhibits a gap for volume fractions ranging from 0.25 to 0.43 (—0.90 < t < 0.00). The maximum Aui/uio = 6.6% 
is reached approximately for a volume fraction of 31% (t = —0.60). The midgap frequency uio ~ 0.773. The F-RD 
structure has thus a gap of comparable width as the PBG in the P structure, but the midgap frequency is much 
higher. For the inverse F-RD structure, which exists for 0.31 < <fi s < 0.78 (—1.07 < t < 1.10), we also find a gap 
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with an even higher midgap frequency, namely ujq = 0.857. The maximum width of the bandgap is however smaller 
(Aoj/ujo — 4.2%) than the one of the direct F-RD structure. In contrast to what is mentioned in Ref. j^, for this 
example it seems that the increase in symmetry of the dielectric structure does not destroy the PBG property but 
allows structures with a higher genus to have a PBG. 

Finally we study three level surfaces belonging to the cubic symmet ry g roup laid: The I2-Y** (called double gyroid 
by Wohlgemuth et al 32}), C(i2-Y**) (belonging to the G' family |32j ) and L (Lidinoid) [^3 surfaces with genus 
g — 9, g = 25 and g = 33, respectively. The group laid contains all the symmetries of the group /4i32, the space 
group of the G surface, as well as the inversion symmetry [42|. For the I2-Y** surface a family of connected level 
surfaces exists for —5.00 < t < 2.98, corresponding to 4> s values between 0.03 and 0.91. Band gap calculations for 
the direct as well as the inverse structure show no PBG, in accordance with the result reported in Ref. [23] for the 
double gyroid structure. This result could suggest that increasing the symmetry of the dielectric structure destroys 
the appearance of a PBG (see Ref. H^)- The C(Ia-Y**) level surface family exists from a volume fraction 4> s of 0.09 
to 0.72 corresponding to t values between t\ = —2.05 and t 2 = 0.99. The band gap maps for the direct and inverse 
C(i2-Y**) structures are shown in FigQ] The results show that the C(h-Y**) structure exhibits a gap for volume 
fractions from 0.10 to 0.56 (—2.00 < t < 0.25). The maximum Auj/uiq = 19.6% is reached approximately for a volume 
fraction of 25% (t — —1.25). The midgap frequency ojq = 0.796. This structure has two attractive features for the 
production of a real PBG crystal: A large gap and a high midgap frequency. This example also shows that increasing 
the symmetry of the dielectric structure does not necessarily destroy the PBG property but allows structures with a 
higher genus to have a PBG. This property was also seen for the direct and inverse F-RD level surface structures. As 
can be seen from FigQ] the inverse C(l2-Y**) structure has a much smaller PBG than the direct C(±2-Y**) structure, 
namely Au/w = 4.7%, and midgap frequency oj = 0.675 for a volume fraction of 39% (t = 0.50). The L level 
surfaces pinch-off at t\ = —0.35 and t 2 = 0.30 which corresponds to values of <fi s from 0.25 to 0.50. The inverse L 
level structures exist for —0.33 < t < 0.30, corresponding to volume fractions 4> s ranging between 50% and 74%. The 
direct as well as the inverse Lidinoid structures do not show a PBG, which could be due to the relatively high volume 
fractions of the high dieletric material. 

IV. TUBULAR STRUCTURES: RESULTS 

For some t < t\ or t > t 2 the P, D and G level surfaces form simple cubic, diamond and gyroid lattices of disconnected 
units With increasing \t\ the minimal surfaces first shrink and look like tubes about the skeletal graphs. 

Then the tubes shrink in the middle while the ends thicken, until the surface pinches-off. The pinching off behavior 
can be controlled by considering the skeletal graphs of the minimal surfaces using level set methods |37j . This results 
in less shrinking of the tubes in the middle and less thickening of the tubes at the ends so that more cylindrical 
tubular structures can be constructed. In Table VI we give the equations for the direct and inverse dielectric tubular 
structures together with 3D renderings of 2x2x2 unit cubes of the structures for various values of t. Note that the 
equations for the direct and inverse structures are the same. The structures become different through Eq.(2) and 
Eq.(3) for the direct and inverse structure, repectively. The interval of t values is chosen such that the endpoints 
approximate the t values for pinch-off of the surface and such that the range of t values for which the structures show 
a PBG is covered. As can be seen from Tables VI the length of these intervals can be quite different for the various 
tubular structures. The corresponding volume fractions cf> s are also given. Note that tubes of dielectric materials are 
in fact only seen for the inverse dielectric structures. 

We study the band gap maps of direct and inverse tubular periodic dielectric structures generated from the skeletal 
graphs of the P, D, and G surfaces [37J. The P structure has sixfold and the D structure has fourfold junctions, 
whereas the G tubular structure consists of triple junctions. Thus the tubular structures have the same topology as 
their corresponding level surface structures. Hence g = 3, g — 5 and g = 9 for the P, G and D tubular structures, 
respectively. 

Our calculations show that the direct tubular P structure, which exists for volume fractions <j) s between 0.29 and 
0.98 (—19.50 < t < 0.09), has a band gap map that is similar to the one of the P level surface structure. As can be 
seen from Fig|S]the tubular P structure exhibits a gap for volume fractions 4> s in the range 19-36%, corresponding to 
t values between -25.00 and -16.50. Note that for —25.00 < t < —19.50 the structure is disconnected and consists of 
"spheres" on a simple cubic lattice. This also happens for the P level surface structure (see section 3). The maximum 
Aw/wo = 4.6% is reached approximately for a volume fraction of 27% (t — —20.00). The corresponding midgap 
frequency u>q — 0.427. The inverse tubular P structure (0.02 < 4> s < 0.72, —19.5 < t < 0.09), which, apart from the 
joints resembles the circular rod scaffold structure in the simple cubic lattice nas no PBG for a dielectric contrast 
of 11.9. Since it has been shown that the square and circular rod structures in the simple cubic lattice have a small 
PBG for a dielectric contrast of 13 [4q|. our results lead to the conclusion that small changes to the structure (lower 
dielectric contrast and slight shape differences) can make these small gaps disappear. 
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The direct (inverse) tubular G structure pinches off at t\ = — 27. 60(— 27.60) and t% = 0.59(0.64), corresponding 
to volume fractions between 0.03 (0.01) and 0.99 (0.97). FigEI shows the photonic band gap maps of the direct and 
inverse tubular G structures. The band gap maps are similar to the band gap map of the level surface structure 
(see also Fig^). For the direct (inverse) tubular G structure the maximum Aw/wo = 20.4%(19.9%) is reached for a 
volume fraction of <f) s = 17%(17%), which corresponds to t = -24.00(-4.00). The midgap frequency is 0.514 (0.502). 

The band gap map of the tubular D structure is split in two, as can be seen from Figd For <f) s < 0.15 (t < —18.90) 
the structure consists of disconnected dielectric spherical units on a diamond lattice (see also Table VI). This structure 
is known to have a PBG |8j. This diamond lattice structure leads to the upper left part of the band gap map. The 
right lower part of the band gap map comes from the connected tubular structure, which exists for —18.90 < t < 0.88 
corresponding to 4> s values between 0.15 and 0.88. The midgap frequencies of the diamond lattice structure are 
almost twice as large as the midgap frequencies of the connected tubular D structure, which causes the splitting in 
the bandgap map. This does not occur in the case of the P level surface structure although also for this dielectric 
structure the band gap map at the lowest volume fractions <j) s is formed by disconnected "spheres" on a simple cubic 
lattice (see Fignj. The direct tubular D structure has a maximal band gap Auj/ujq — 19.2% with midgap frequency 
ujq = 0.546 for a volume fraction of (j> s — 25%, which corresponds to t = —17.00. Note that the band gap map 
of the corresponding D level surface structure (see Fig^l is not split in two. This is due to the fact that only the 
connected structures have a PBG for certain volume fractions and that at pinch-off, the D level surface structure does 
not become a diamond lattice structure. 

The inverse tubular D structure has a band gap map that is similar to the band gap map of the corresponding 
level surface structure, as can be seen from FigQ] and Fig[7| For the inverse D tubular structure the maximum 
Aw/wo = 21.4% is reached for a volume fraction of cj> s = 22%, which corresponds to t — —4.00. The midgap frequency 
is 0.569. The results are summarized in Table VII. 



V. CONCLUSIONS 

For a fixed dielectric contrast the presence of a photonic band gap is due to a complex interplay of composition 
(volume fraction of structure and background), the symmetry and topology of the periodic dielectric structure. In 
general, within one symmetry group, the cubic structures with a simple topology (low genus) show a band gap and the 
complex structures with a high genus do not. Increasing the symmetry of the dielectric structure does not necessarily 
destroy the photonic band gap property but allows structures with a higher genus to have a photonic band gap. 

For the production of photonic band gap materials the D, G, and C(l2-Y**) structures are the most interesting 
ones. They all have large gaps of approximately 20%. For the low wavelength regime the C(l2-Y**) structure with 
IaZd symmetry might be most favourable since its midgap frequency (uq = 0.796) is a factor of 1.5 higher than the 
one for the (tubular) D and G structures. 

The tubular D and G structures may be easier to produce than their level surface structure counterparts and have 
similar band gap widths and midgap frequencies. 
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TABLE I: Equations of various level surfaces [32I |3^| and 3D renderings of 2x2x2 unit cubes of level surface structures for 
various values of £ and corresponding volume fractions S . M = 2irm/a where M = X,Y,Z, m = x,y,z and a denotes the 
length of the crystallographic cell. 

C( ± Y) -2cosXcosYcosZ + sin2XsinY + sin X sin 2Z + sin 2Y sin .Z = £ 




t = —1.6 £ = -1.0 £=-0.5 t = 0.0 £ = 0.5 t = 1.0 t = 1.6 

S = 0.08 (j> s = 0.25 <j) s = 0.38 S = 0.50 S = 0.62 S = 0.75 S = 0.92 



± Y 2 cos X cos F cos Z + sin2XsinF + sinXsin2Z + sin2YsinZ = £ 

^Sk *si sffli -^s 

&5tP : .v*P -Cs-p % 

£ = -0.4 £ = -0.3 £ = 0.0 £ = 0.3 £ = 0.4 

<t> s = 0.34 S = 0.37 S = 0.50 S = 0.63 S = 0.66 

C(S) cos 2X + cos 2Y + cos 2Z + 2[sin 3X sin 2Y cos Z + cos X sin 3Y sin 2Z + sin 2X cos Y sin 3Z] 
+2[sin 2X cos 3 Y sin Z + sin X sin 2 Y cos 3Z + cos 3X sin Y sin 2Z] = £ 




£ = -1.4 £ = -1.0 £ = -0.5 £ = 0.0 £ = 0.5 £ = 1.0 £ = 1.4 

4> s = 0.24 S = 0.30 S = 0.40 S = 0.50 S = 0.60 S = 0.70 S = 0.76 



C(Y) — sin X sin Y sin Z + sin 2X sin Y + sin 2Y sin Z + sin X sin 2Z — cos X cos Y cos Z 
+ sin 2X cos Z + cos X sin 2F + cos Y sin 2Z = £ 




£ = -1.8 £ = -1.5 £ = -1.0 £ = 0.0 £=1.0 £ = 1.5 £ = 1.8 

4> s = 0.09 0s = 0.16 S = 0.28 4> s = 0.50 S = 0.72 S = 0.84 S = 0.91 



Y cos X cos Y cos Z + sin X sin Y sin Z + sin 2X sin Y + sin 2 Y sin Z + sin X sin 2Z 
sin 2X cos Z + cos X sin 2 Y + cos Y sin 2Z = £ 




£ = -0.4 £ = 0.0 £ = 0.4 

4> s = 0.29 0s = 0.50 4> s = 0.71 



G sin X cos Y + sin Y cos Z + cos X sin Z = t 




£ = -1.4 £ = -1.0 £ = -0.5 £ = 0.0 £ = 0.5 £ = 1.0 £ = 1.4 

4> s = 0.02 0s = 0.17 0s = 0.34 0s = 0.50 S = 0.66 S = 0.83 S = 0.98 
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TABLE II: Same as Table 1 
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TABLE HI: Same as Table 1 

F-RD 4 cos X cos Y cos Z — [cos 2X cos 2 Y + cos 2X cos 2Z + cos 2Y cos 2Z] = t 




t = -1.0 t = -0.5 t = 0.0 £ = 0.5 t = 0.7 

4> s = 0.23 4> s = 0.33 4> s = 0.43 4> s = 0.54 <f) a = 0.60 



inv F-RD 4 cos X cos Y cos Z — [cos 2X cos 2 Y + cos 2X cos 2Z + cos 2Y cos 2Z] = £ 




£ = -1.0 t = -0.5 £ = 0.0 £ = 0.5 £=1.0 

</> s = 0.77 S = 0.67 S = 0.57 S = 0.46 <f> a = 0.33 



D sin X sin Y sin Z + sin X cos Y cos Z + cos X sin Y cos Z + cos X cos Y sinZ — t 




£ = -1.0 i = -0.5 £ = 0.0 £ = 0.5 £=1.0 

S = 0.16 S = 0.33 (j) s = 0.50 4> s = 0.67 S = 0.84 



C(D) cos(3X + Y) cos Z - sin(3X - Y) sin Z + cos(X + 3Y) cos Z + sm{X - 3Y) sin Z 
+ cos(X - Y) cos 3Z - sin(X + Y) sin 3Z = £ 




£ = -0.3 £=-0.1 £ = 0.0 £ = 0.1 £ = 0.3 

4> s = 0.34 4> s = 0.40 4> s = 0.46 S = 0.56 <t> a = 0.59 



I 2 _ Y" -2 [sin 2X cos Y sin Z + sin X sin 2Y cos Z + cos X sin Y sin 2Z] 
+ cos 2X cos 2Y + cos 2 Y cos 2Z + cos 2X cos 2Z = £ 




£ = -5.0 £=-4.0 £=-3.0 £ = -2.0 £ = 0.0 £ = 2.0 

4> s = 0.03 (f} s = 0.10 S = 0.18 4> s = 0.26 <j) a = 0.43 S = 0.68 



inv h - Y** -2 [sin 2X cos Y sin Z + sin X sin 2Y cos Z + cos X sin Y sin 2Z] 
+ cos 2X cos 2Y + cos 2Y cos 2Z + cos 2X cos 2Z = £ 




£ = -4.0 £=-3.0 £=-1.0 £ = 0.0 £ = 1.0 £ = 2.0 £ = 2.9 

4> s = 0.90 (j) s = 0.82 4> s = 0.66 S = 0.57 S = 0.46 <f> s = 0.323 <j) a = 0.10 
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TABLE V: Maximum band gap widths Aoj/cjo, corresponding midgap frequency u>o, volume fraction <f) s and parameter t of 
various single level surface structures, categorized by their space group and genus g. A dash indicates that the structure has no 
photonic band gap. a The value found in the literature for the genus for the crystallographic cell of the minimal C(S) surface is 
g = 17 |43]l . 6 The value found in the literature for the genus for the crystallographic cell of the minimal Y surface is a = 17 |43|. 
c Same value as reported by Garstecki et al f44j, but different from the value g = 145 reported by Fisher et al. |43j| for the 
minimal D surface. For an explanation of the deviating values of the genus see text. 



Space group (index) 


Surface 


9 


4>s (%) 


t 


Alu/uj (%) 


LOO 


Po3 (205) 


CPY) 
±Y 


13 
21 


17 
- 


-1.30 
- 


8.1 

- 


0.786 
- 


Ja3 (206) 


C(S) 


65 a 










P4 3 32 (212) 


C(Y) 


13 


22 


-1.25 


6.2 


0.630 




Y 


13 b 










74i32 (214) 


G 


5 


21 


-0.90 


21.7 


0.474 




C(Y**) 


29 


18 


-3.25 


11.9 


0.514 


I43d (220) 


S 


21 


_ 


- 


- 


- 


Pm3m (221) 


P 
I-WP 
inv I-WP 


3 
7 
7 


24 
. 


-0.90 


6.8 


0.442 




C(P) 


9 










Fm3m (225) 


F-RD 


21 


31 


-0.60 


6.6 


0.773 




inv F-RD 


21 


39 


0.80 


4.2 


0.858 


Fd3m (227) 


D 

C(D) 


9 
121 c 


25 


-0.60 


20.3 


0.541 


Ia3d (230) 


I 2 -Y** 
inv I 2 -Y" 


9 
9 












C(I 2 -Y") 


25 


25 


-1.25 


19.6 


0.796 




inv C(I 2 -Y**) 


25 


39 


0.50 


4.7 


0.675 




L 


33 












inv L 


33 











[45] H.S. Soziier, and J.W. Haus, J. Opt. Soc. Am. B 10, 296 (1993). 
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TABLE VI: Equations and 3D renderings of 2x2x2 unit cubes of various tubular structures for various values of t and corre- 
sponding volume fractions <fr s . M = 2nm/a and M' — M — 7r/4, where M = X, Y, Z, m — x,y, z and a denotes the length of 
the crystallographic cell. 

tubular P 10[cos X + cos Y + cos Z] — 5.1[cosX cos Y + cos Y cos Z + cos Z cosX] — 14.6 = t 




t = 0.0 
= 0.97 



inv tubular P 10[cosX + cos Y + cos Z] — 5.1[cosXcos Y + cos Y cos Z + cos Z cos X] — 14.6 = t 




t = 0.0 
(j> s = 0.03 



tubular D 10 [sin X' sin Y' sin Z' + sin X' cos Y' cos Z' + cos X' sin Y' cos Z' + cos X' cos Y' sin Z'] 
-0.7[cos4X + cos4y + cos4Z] - 11 = t 




inv tubular D 10 [sin X' sin V" sin Z' + sin X' cos Y' cos Z' + cos X' sin F' cos Z' + cos X' cos F' sin Z'] 
-0.7[cos4X + cos4y + cos4Z] - 11 =t 




t = -20.0 t=-15.0 t=-10.0 t=-5.0 t = 0.0 t = 0.5 

S = 0.99 S = 0.65 S = 0.45 <f> a = 0.26 S = 0.06 <j> s = 0.03 



tubular G 10[cos X sin F + cos Y sin Z + cos Z sin X] - 0.5[cos 2X cos 2Y + cos2Fcos2Z + cos2Zcos2X] - 14 = t 




t = -25.0 t=-20.0 t=-15.0 t = -10.0 t = 0.0 t = 0.5 

S = 0.13 <j) s = 0.31 S = 0.47 (j)s = 0.63 <f> s = 0.97 <^) s = 0.99 



inv tubular G 10[cos X sin Y + cos Y sin Z + cos Z sin X] - 0.5[cos 2X cos 2Y + cos 2V cos 2Z + cos 2Z cos 2X] - 14 = t 




t = -20.0 t=-15.0 t=-10.0 £=-5.0 t = 0.0 t = 0.5 

S = 0.69 3 = 0.53 4> s = 0.37 S = 0.21 S = 0.03 4> 3 = 0.01 
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TABLE VII: Maximum band gap widths Aw/wo, corresponding midgap frequency uio, volume fraction 4> s and parameter t of 
various tubular structures, categorized by their genus g. A dash indicates that the structure has no photonic band gap. 



Tubular structure 


9 


S (%) 


t 


Aw/w„ (%) 




P 


3 


27 


-16.50 


4.6 


0.427 


inv P 


3 










G 


5 


17 


-24.00 


20.4 


0.514 


inv G 


5 


17 


-4.00 


19.9 


0.502 


D 


9 


25 


-17.00 


19.2 


0.546 


inv D 


9 


22 


-4.00 


21.4 


0.569 



